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0\ ■ The density of states of a Bose-condensed gas confined in a harmonic trap 
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is investigated. The predictions of Bogoliubov theory are compared with 



. the ones of Hartree-Fock theory and of the hydrodynamic model. We show 

C 

that the Hartree-Fock scheme provides an excellent description of the exci- 



tation spectrum in a wide range of energy, revealing a major role played by 

single-particle excitations in these confined systems. The crossover from the 

hydrodynamic regime, holding at low energies, to the independent particle 

regime is explicitly explored by studying the frequency of the surface mode 

as a function of their angular momentum. The applicability of the semiclas- 

sical approximation for the excited states is also discussed. We show that the 

semiclassical approach provides simple and accurate formulae for the density 

of states and the quantum depletion of the condensate. 
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I. INTRODUCTION 



The collective modes of a Bose condensed gas confined by an external potential have 
been the object of extensive work in the last months. The successful agreement between 
experimental results and theoretical predictions for the collective frequencies at 
low temperature has stimulated an intensive research activity. Though only the modes 
with low multipolarity and frequency have been detected in experiments, the excitations 
at higher energy and angular momentum are also very important because they determine 
the statistical behavior of the system, including thermodynamics, transport phenomena and 
superfluid effects. 

The excited states at high energy are expected to have single-particle nature. However, 
the transition from collective (phonon-like) to single-particle excitations in a inhomogeneous 
system can differ significantly from the case of a uniform Bose gas. In fact, the presence of 
a surface allows for the occurrence of single-particle states even at low energy (lower than 
the chemical potential). These states, of low energy but high multipolarity, are localized 
near the surface, where the condensate density becomes small. This behavior represents a 
peculiar and interesting feature of these confined systems; in a uniform Bose gas, in fact, 
only phonons are present at low energy. In a recent paper |9j we have already pointed out 
the effects of these single-particle states on the thermodynamic properties of the trapped 
gases. 

In the present work we solve the equations for the excited states of a weakly interacting 
gas in a spherical trap at zero temperature within Bogoliubov theory. The main purpose 
is to investigate the collective (phonon-like) and single-particle character of the elementary 
excitations. This is accomplished by calculating key quantities, such as the density of states, 
the frequency of the surface modes and the quantum depletion of the condensate, and by 
comparing the predictions of Bogoliubov theory with the ones of different approximations, 
like Hartree-Fock theory and the hydrodynamic model. Finally we check the accuracy of 
the semiclassical approximation and we show that it provides simple and useful formulae for 
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both the density of states and the quantum depletion of the condensate. 



II. BOGOLIUBOV THEORY 

The elementary excitations of a degenerate Bose gas are associated with the fluctuations 
of the condensate. At low temperature they are described by the time dependent Gross- 
Pitaevskii (GP) equation for the order parameter [ID|: 



i^*(r, t) = [-~ + V ext (r) + g | *(r, t) | 2 ) *(r, t) , (1) 

where Jc?r|\l/| 2 = iV is the number of atoms in the condensate. At zero temperature N 
coincides with the total number of atoms, except for a very small difference 5N <C iV due 
to the quantum depletion of the condensate. The coupling constant g is proportional to the 
s-wave scattering length a through g = 4irh 2 a/m. In the present work we will discuss the 
case of positive scattering length, a > 0, as in the experiments with Rubidium and Sodium, 
but the same formalism can be also applied to systems with negative scattering length. The 
trap is included through V ex t, which is chosen here in the form of an isotropic harmonic 
potential: V ex t(r) = {l/2)muj 2 HO r 2 . The harmonic trap provides also a typical length scale 
for the system, clho — (^/mw/zo) 1 ^. Actually, the experimental traps have cylindrical 
symmetry, with different radial and axial frequencies, but the choice of a spherical trap, as 
we will discuss later, is not expected to affect the main conclusions of the present work, 
while reducing greatly the numerical effort. 

The normal modes of the condensate can be found by linearizing the GP equations, i.e., 
looking for solutions of the form 

#(r, t) = e~ ilMt/h [tt„(r) + u(r) e - <w * + v*(r)e iut ] (2) 

where \i is the chemical potential and functions u and v are the "particle" and "hole" com- 
ponents characterizing the Bogoliubov transformations. After inserting in (JI|) and retaining 
terms up to first order in u and v, one finds three equations. The first one is the nonlinear 
equation for the order parameter of the ground state flOp , 



[/To + ^(r)]* (r) = ^*o(r), (3) 

where Hq = — (h 2 /2m)'V 2 + V ex t(v); while u(r) and t>(r) obey the following coupled equation: 

hwu(r) = [H -/2 + 2g^ 2 }u(r) + g* 2 v(r) (4) 
-huuv{r) = [H -fi + 2gVl]v{r) + gV 2 u(r) . (5) 



Numerical solutions of these equations have been recently found by different authors 
In the present work, we use them to calculate the density of states, the frequency of the 
surface modes, the quantum depletion of the condensate, in order to clarify the different role 
played by excitations having collective and single-particle character. 

When the dimensionless parameter Na/a HO is large, the kinetic energy term in the 
ground state equation (|3|) becomes negligible with respect to the mean field term and one 
gets the Thomas-Fermi approximation: 

*7{r)=( ^ F ~ Vext{r) \ 12 (6) 
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In the same limit the equations of motion coincide with the equations of the hydro- 
dynamics (HD) of superfluids P, |TT| , |i"2|] . In the spherical case their eigenfrequencies take the 
analytic form 0] 

u(n,£) = to H o(2n 2 + 2n£ + 3n + i) 1/2 , (8) 

where I and n are the angular momentum quantum number and the number of nodes in 
the radial solution, respectively. The deviations from the predictions of the noninteracting 
harmonic oscillator (HO) model, 

uo(n,£) = LU HO (2n + £) } (9) 

point out the effects of two-body interactions. These are particularly important for the so 
called "surface" modes (n = 0), where the HO prediction u = £u>ho is significantly lowered 



to the hydrodynamic value uj = y/I ujho- in general the HD prediction turns out to be 
very accurate for the low energy excitations of the system, while the ideal gas prediction is 
expected to be valid in the opposite case of high excitation energies. The exact solutions of 
the equations (§-[1) provide the correct interpolation between the two limiting regimes. 

A typical spectrum obtained from equations (§-0) is given in the upper part of Fig. [I] for 
a gas of N = 10000 atoms of Rubidium (scattering length a = llOao, where ao is the Bohr 
radius). For the spherical trap we have chosen the frequency ujho = 2%uho = 27rl87Hz, 
which is the average to ho — (^x^y^z) 1 ^ 3 of the axial and radial frequencies of Ref. Q. It 
corresponds to the oscillator length ano — 0.791 x 10 -4 cm. Energy is given in units fuvno 
and the chemical potential is 8.41 in these units. The vertical bars have length (21 + 1), 
so that the angular momentum of each state can be inferred from the figure. One clearly 
sees that, at energy much larger than the chemical potential, the excited states tend to be 
grouped into levels Huoho apart, as in the noninteracting HO model. Conversely, the energy 
of the lowest modes is close to the prediction of the HD equations ||. For instance, the 
lowest £ = 2 and £ = modes differ by less than 3% from the hydrodynamic values a/2 and 
y/E, respectively. 



III. DENSITY OF STATES 

Once the spectrum of excited states is calculated, one can count the number of states 
below a given energy e: 

W(e) = X;(2*+l). (10) 

The density of states, g(e) is simply the derivative d%l(e)/de. Note that the quantity 91(e) 
is well defined even for a discretized spectrum, while the density of states implies averaging 
the number of states within small but finite energy intervals. 

In Fig. |2| we show the quantity 91(6) obtained by counting the levels in the spectrum of 
Fig. [J]. For comparison the results of the noninteracting harmonic oscillator (^) and of the 
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hydrodynamic model (j8|) are also shown. The effects of the repulsive interatomic forces are 
clearly responsible for an enhancement of the density of states with respect to the ideal gas. 
However 9?(e) remains well below the HD approximation, the latter being soon inadequate 
as e increases. Indeed hydrodynamic theory accounts for collective phenomena and provides 
an excellent description of the low-lying elementary excitations of the system |3|,|3|], but 
completely ignores single-particle effects. This is exactly the opposite of what Hartree-Fock 
(HF) theory does. For this reason it is interesting to compare the results of Bogoliubov 
theory with the predictions of HF theory in which one determines the eigenstates of the 
single-particle Hamiltonian |T3|JT^ 



Hhf = -(h 2 /2m)\/ 2 + V ext {v) + 2<?^(r) . (11) 

The lowest eigenstates of the HF Hamiltonian are expected to be localized near the surface 
of the condensate. To understand better this point, let us take the large N limit. In this 
case one can use the Thomas-Fermi approximation (Q) for the ground state density. The 
HF Hamiltonian then takes the simple form 

Hhf = -{h 2 /2m)V 2 + ^mu 2 HO \r 2 - R 2 \ , (12) 

where R = [2/i TF /(mw^ )] 1/ ' 2 is the classical radius of the condensate. The HF potential 
has a pronounced minimum at R. This potential well near the boundary persists in the HF 
Hamiltonian even for smaller values of N. 

It is worth stressing that, in general, Hartree-Fock theory is expected to be correct for 
energies larger than the chemical potential. For these trapped bosons, however, it accounts 
also for the low energy excitations close to the boundary, where the density of the condensate 
is small. This can be seen in Fig. [I], where the Hartree-Fock and Bogoliubov spectra are 
compared. One notes that the two spectra are rather similar even below the chemical 
potential (/z = 8.41 in this case) except for the excitations having lowest angular momenta, 
i.e., shortest bars in the figure. Those collective excitations are instead correctly reproduced 
by the hydrodynamic model. The corresponding Hartree-Fock prediction for the quantity 
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91(e) is also given in Fig. [| as a solid line. The agreement with the results of equations (f§-[5D 
is remarkable in the whole range of energy. The figure shows the case of 10000 atoms but a 
similar agreement is found for all values of N relevant for the experiments (N ~ 10 3 -j- 10 7 ). 
The above behavior represents a major difference with respect to the case of a homogeneous 
Bose gas where phonons have a crucial effect on the density of states. 



IV. SURFACE EXCITATIONS 

In order to better understand the transition from the collective to the single particle 
regime, we have explored in details the evolution of the excitation energy for the surface 
modes (n = 0) as a function of their angular momentum £, as predicted by the solution of 
the equations As already mentioned, the effects of the interaction are particularly 

important for such modes and are responsible for a significant lowering of their frequency. 
In Fig. we plot the quantity uii/i, in units ujhoi f° r different numbers of atoms in the same 
trap. This ratio has an important physical meaning because, according to Landau's criterium 
for superfluidity, it provides the rotational frequency at which the £-th surface excitation 
becomes unstable. The hydrodynamic prediction (H) is shown as a dashed line, the curve 
going asymptotically to zero for large £. The figure shows that the Bogoliubov states first 
follow the HD curve, but, rather soon, they deviate from it, approaching asymptotically the 
non interacting value u>e/£ = who- The deviation from HD takes place at larger values of I if 
N is increased, revealing that the HD approximation becomes applicable to a larger number 
of states in this limit. 

A simple estimate of the typical value of i at which the HD picture starts failing, one can 
take |L5| £ c ~ Rp c where R is the radius of the condensate, proportional to A 1//5 , and p c is of 



the order of the inverse of the surface thickness d = [a% / (2-R)] 1 ^ 3 |I6|-|I8||. For larger values 
of £, the wavelength of the excitations becomes shorter than d and one explores microscopic 
details of the boundary which can not be described by the Thomas-Fermi approximation 
(0) and by the HD equations. This yields l c oc (R/a HO )^ 3 oc A 4//15 , corresponding to an 
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excitation energy hw oc iV 2 / 15 , smaller than the chemical potential, which instead behaves 
as N 2 / 5 . This explains why the crossover from the HD to the single particle regime takes 
place at energies smaller than \i. 

For each value of N, the curves in Fig. |^ exhibit minima and one can define a critical 
frequency as Q c = mm(ut/£). For rotational frequencies larger than Q c the surface excita- 
tions become unstable. It is interesting to compare this value with the critical frequency 
needed to generate a vortex [BJ] . This is done in Fig. [| where we compare the two critical 
frequencies as a function of N. We find that the lowest instability is always associated with 
the creation of a vortex. Note however that in order to generate a vortex one needs to 
transfer to the system a huge angular momentum (equal to Nh) , which is much higher than 
the value £ required to create a surface excitation. 



V. SEMICLASSICAL APPROXIMATION AND SCALING BEHAVIOR 



A good approximation for the density of states can be obtained by solving equations 
in the semiclassical approximation 



, pT| . In this approximation, which is expected 
to hold for excitation energies much larger than the oscillator energy Hujho, the quantity 
%l(e) is a continuous function of e defined by 

drdp 



de 



6(e'-e(r,p)) 



(13) 



where 



e(r,p) 



p 

2m 



+ V ext (r) - fi + 2gn (r) 



9 X( r ) 



1/2 



(14) 



corresponds to the semiclassical dispersion law. Here the quantity no(r) = ^(r) is the 
condensate density. In Fig. |^a we compare the semiclassical result for 9T(e)/iV (solid line) 
with the one obtained from equations (§-^) (squares) for 10 4 atoms of Rubidium in the same 
trap of Fig. ||. Here the energy is given in units ksT c = huiHoiN/Cfi)] 1 / 3 , which is the 
critical temperature for an ideal Bose gas in a harmonic trap; the value for 10 4 atoms is 



S 



ksT c = 20.26Hujho as shown also in Fig. |I[ The accuracy of the semiclassical approximation 
turns out to be very high also for relatively low values of e. 

The use of the semiclassical approximation allows one to carry out the analysis of the 
density of states in a systematic way and to exploit the dependence on the relevant param- 
eters of the system. In fact, when the number of atoms in the condensate is large enough to 
make the Thomas-Fermi approximation (P) accurate, the statistical properties of the system 
can be expressed in terms of a single scaling parameter 77 given by the ratio |22| 



k B T c \ ano J 

between the chemical potential (^), calculated at zero temperature in the Thomas- Fermi 
approximation, and the critical temperature /cbT c . The ratio rj depends on the defor- 
mation of the trap only through the geometrical average of the oscillator frequencies 
^ho = (wjWyW^ 1 / 3 . The parameters used in the calculation of Fig. [|a correspond to a 
spherical trap with rj = 0.407. As pointed out in Ref. p2| , quite different experimental 
conditions (shape of the trap, value of N, etc.) can correspond to very similar values of 77. 
In terms of the scaling parameter 77 and the dimensionless energy e = e/(A%T c ), the number 
of states 9t(e) predicted by the Bogoliubov semiclassical theory becomes 



N Jo vrC(3) 



d e — rrzz / dx y/l — X 
Jo 



X 



e 77— — . = + (0 V « + We 

y/x 2 + (e') 2 /ri 2 



(16) 



This result has been obtained by using the Thomas-Fermi approximation (|^) for the con- 
densate density in ([T3| - |T4]) ; this allows one to split the space integral into an inside region 
(first term in the square bracket) and an outside region (second term). It is worth stressing 
that equation ([16]), which are expected to hold in the scaling regime Na/ano ^ 1> provide 
a very good estimate of the semiclassical expression (0) even for relatively small N. For 
instance, the two predictions are indistinguishable in Fig [5]a, being represented by the same 
solid line. In Fig. [5|b we show the predictions for the density of states given by the semi- 
classical approximation fll6|) for three different values of 77. The parameters of the recent 
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experiments at Jila |l] and MIT |J, using very different traps, correspond to 77 ranging from 
0.39 to 0.45. 

Expression flTED can be also expanded at low energy, e <C ksT c , still compatible with the 
assumption e ^> fkono- One finds the law 9t(e)/iV oc e 5 / 2 . This differs from the usual e 3 law 
typical of the phonon regime, revealing the different behavior exhibited by these systems 
with respect to the homogeneous Bose gas. 



VI. QUANTUM DEPLETION 

In the last part of the paper we calculate the quantum depletion of the condensate which, 
according to Bogoliubov theory, is given by 

3 

The "hole" components Vj can be obtained by solving equations (|]-[5|). In semiclassical 
approximation || one replaces the sum over all the discrete states with the integral over p 
of the function 

V2{P) r) = 2^Trj (fk + Vextiv) ~ ^ + 29710 ^ + C(P ' r) ) ' (18) 



where e(p, r) is the single particle energy fll4D - In a uniform gas this expression yields 
the most famous result SN/N = (8/3)(n a 3 /n) 1 ^ 2 . In the trapped gas and in the limit 
Na/dHo ^ 1) where TF approximation holds, the semiclassical approximation provides the 
simple analytic law M 



5N 7] 
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0.098 T] 6 , (19) 



N 6V2C(3) 

with i] given in Eq. (|15|). Since the available experiment correspond to 77 ~ 0.4, the quantum 
depletion turns out to be less than 1%, as already pointed out in Refs. HJ^]. In Fig. |] we 
show the quantum depletion for 10000 and 50000 atoms of Rubidium obtained by summing 
over the Bogoliubov spectrum up to a given energy e (solid lines). We compare it with the 
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prediction of the semiclassical expression fll8|) (dashed lines), while the arrows indicate the 
asymptotic values (fTDD, holding in the scaling limit. An important result emerging from the 
numerical calculation is the very slow convergency of the sum fll?p . This is not a surprise, 
since also in a homogeneous gas the convergency is slow due to the 1/p 4 tail in the momentum 
distribution at high momenta and one has to go up to e = lOO/i in order to saturate 90% of 
the sum fll7|). 

The agreement between the quantum depletion obtained from the discretized sum fll7|) 
over the Bogoliubov states and from the semiclassical approximation is satisfying and 
was not obvious a priori. Fig. || shows a discrepancy of the order of 5% between the two 
predictions for iV = 10000, while for larger N the two curves tend to coincide. It is worth 
noticing that the two solid lines in Fig. |^ requires the summation of J dr\vj(r)\ 2 over up to 
15000 different values of (n, £) in the Bogoliubov spectrum; the calculation is then much 
heavier than the semiclassical one. The good accuracy of the semiclassical approach makes 
it useful in practical situations. This is especially true for the simple formula (|1^) which 
includes the case of anisotropic traps through the averaged frequency to ho — {^x^y^z) 
entering the Thomas-Fermi chemical potential /i TF and, hence, the scaling parameter rj. 

Finally, the large N semiclassical formula ( |T9"D shows the rather strong dependence of 
the depletion on the scattering length parameter a. If the magnetic tuning of the scattering 
length will become available, it will be possible in the future to increase significantly the 
value of 7] and consequently explore Bose gases where the quantum depletion is much larger. 

VII. SUMMARY 

We have investigated the elementary excitations of a dilute Bose gas in harmonic trap 
by solving the equations of Bogoliubov theory. Differently from the case of a uniform gas, 
where phonons dominate the system at energies of the order of, or lower than the chemical 
potential, the spectrum of the trapped gas shows an important particle-like behavior even at 
low energy. This fact has been here explored in detail. We have compared the results of Bo- 
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goliubov theory for the density of states with the ones of Hartree-Fock theory finding a very 
good agreement on a wide range of energy. We have studied the behavior of surface modes, 
emphasizing the crossover from the low energy regime, well described by the hydrodynamic 
model, to the single-particle regime. This crossover provides also a critical frequency as- 
sociated with a rotational instability and we have compared this frequency with the one 
needed to create a quantized vortex. Another important result emerging from our analysis 
is the high accuracy exhibited by the semiclassical approximation for the excited states. 
Finally, we have calculated the quantum depletion of the condensate by summing the "hole" 
component fdr\v(r)\ 2 over all the states in the excitation spectrum of Bogoliubov theory. 
The convergency of the sum turns out to be very slow, as expected by the analogy with the 
case of the uniform gas. Again we find an excellent agreement with the predictions of the 
semiclassical approximation. In the limit Na/ano ^ 1, the latter provides the simple and 
useful formula 5N/N = 0.098 77 s , in terms of the scaling parameter 77 = l.57(N 1 / 6 a/aHo) 2 ^- 
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FIGURES 

FIG. 1. Excitation spectrum of 10000 atoms of 87 Rb in a spherical trap with clho = 0.791 x 10~ 4 
cm. The vertical bars have length (21 + 1). The upper spectrum corresponds to the numerical 
solution of equations (|4||5|); the lower one is the spectrum of the Hartree-Fock Hamiltonian (11). 



Two energy scales (in units of friono) are a l so shown in the figure: the chemical potential \x = 8.41, 
fixed by the solution of the Gross-Pitaevskii equation (||), and the critical temperature for a non- 
interacting gas in the same trap, ksT c = 20.26. 

FIG. 2. Number of states 91(e) vs. energy. The Bogoliubov (points) and Hartree-Fock (solid 
line) predictions, obtained by counting the states in Fig. [I], are compared with the ones of the 
noninteracting harmonic oscillator (dashed line) and of hydro dynamic equations (dot-dashed line). 

FIG. 3. Frequency (in units of ujho) of the n = excited states as a function of their an- 
gular momentum i for N atoms of 87 Rb in a spherical trap with clho = 0.791 x 10~ 4 cm. The 
hydrodynamic prediction is shown dashed line. 

FIG. 4. Critical rotational frequency (in units of coho) f° r producing of a quantized vortex 
(solid line) or surface states (dashed line) as a function of the number of Rubidium atoms in the 
spherical trap. 

FIG. 5. Ratio 9?(e)/iV vs. e, in units /c#T c . In part a, the squares correspond to counting the 
states in the Bogoliubov spectrum of Fig. |], while the solid line is the corresponding semiclassical 
approximation (|l~3|). The latter is indistinguishable from the formula (|i~6|), valid in the scaling 
regime Na/auo S> 1- In part b, the semiclassical prediction ( |l6|) is given for different values of the 
scaling parameter r\. 
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FIG. 6. Quantum depletion for 10000 (two lower curves) and 50000 (two upper curves) atoms 
of Rb in a spherical trap with ano = 0.791 x 10~ 4 cm. The depletion is plotted as a function of 
the maximum energy considered in the sum (|l7|). Solid lines: \v\ 2 from the solution of equations 
|]5|); dashed lines: from the semiclassical approximation Qi~8|). Arrows: asymptotic scaling values 
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